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摘 要 :城市 土壤 有 机 碳 (SOC ) 分 布 受 城 市 建设 工业 发 展 等 人 为 因素 的 影响 表现 出 明显 的 空间 差 
异 。 为 揭示 石嘴山 市 SOC 受 城 市 化 .工业 化 等 人 类 活动 的 影响 ,分 别 利 用 普通 克 里 格 法 (OK)、 多 元 
线性 回归 克 里 格 法 (RK) 3E Ek CZ 3& (RS) FE gk — 8 38 Jo C [el 98 v6 8 A (RGWRK) 83 6 "E 
山 市 SOC 空间 分 布 。 结 果 表 明 : 石 嘴 山 市 SOC 含量 在 1.31 ~ 6692 g-kg :之 间 变 化 ,其 平均 值 为 
17.61 gkg!'。 石 嘴 山 市 不 同 功能 区 SOC 含量 存 在 显著 差异 (p < 0.05) ,具体 表现 为 工业 区 > 医疗 区 > 


商业 区 > 道路 > 住宅 区 > 公园 > 农田 > 科教 区 ;SOC 含 


量变 异 系数 为 66.27% , 呈 中 等 程度 变异 ;其 最 佳 


拟 合 模型 为 高 斯 模型 ,CV(Co+C) 为 0.02, 属 于 强 空 间 自 相关 。SOC 与 遥感 影像 波段 DN 值 的 差 值 
(B1-B7、B3-B7、B4-B7) 和 和 地形 因子 (高 程 .坡度 、 起 伏 度 ) 之 间 存 在 着 极 显著 的 相关 性 (p < 0.01); 
通过 对 4 种 方法 的 结果 进行 对 比 可 知 以 各 波段 DN 值 差 值 与 地 形 因 子 为 输入 量 , 利 用 RGWRK 预测 
的 SOC 精度 最 高 , 相 比 OK 精度 提高 了 10.05% , 相 比 RK 精度 提高 了 8.79% , 相 比 RS 精度 提高 了 
8.92%; 研 究 区 SOC 含量 呈 北 高 南 低 的 趋势 ,工业 区 SOC 含量 是 郊区 农田 SOC 含 量 的 1.92 们 ,表明 


城区 SOC 含量 有 富 集 的 趋势 。 


关键 词 : 城市 土壤 有 机 碳 ; 空间 预测 ; 遥感 反 演 ; 地 理 加 权 回 归 克 里 格 


文章 编号 : 


土壤 是 参与 地 球 与 大 气 碳 循环 的 陆地 生态 系 
统 的 重要 组 成 部 分 ,并 对 全 球 环境 的 稳定 起 着 至 关 


提供 了 科学 依据 。 
遥感 技术 是 可 以 在 短 时 间 内 对 地 面 同一 地 区 


重要 的 作用 "。 其 中 城市 土壤 是 一 个 巨大 的 碳 库 ， 
其 碳 循环 是 全 球 碳 循环 的 重要 组 成 部 分 。 随 着 工 
业 发 展 和 城市 扩张 , 越 来 越 多 的 自然 和 农业 用 地 正 
逐步 变 为 城市 ,这 将 对 全 球 气候 的 变化 和 陆地 生态 
系统 的 稳定 产生 深远 的 影响 ”。 土 壤 有 机 碳 (SOC) 
是 土壤 质量 的 核心 ,是 土壤 碳 与 大 气 C0: 快 速 交换 
的 主要 形式 ”。 在 城市 环境 背景 下 ,由 于 土地 利用 
方式 的 改变 和 人 类 频繁 的 活动 ,使 得 城市 SOC 含量 
存在 较 大 的 空间 变异 性 。 因 此 ,城市 化 进程 中 SOC 
含量 和 动态 变化 趋势 的 高 效 准 确 分 析 和 预测 ,为 城 
市 土壤 的 管理 \ 利 用 ,公共 卫生 和 城市 环境 的 维护 
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实现 反复 观测 的 一 种 新 兴 综 合 性 勘探 技术 ,同时 可 
以 获得 大 量 的 和 遥感 影像 数据 ,这 些 遥 感 影像 记录 了 
土壤 的 某 些 属性 ,在 土地 评价 中 应 用 广泛 外 。 目 前 
应 用 于 多 光谱 遥感 影像 SOC 的 预测 方法 主要 有 多 
元 线性 回归 外 最 小 二 乘 回归 等 ,这 些 方法 在 特定 
的 研究 区 均 取 得 了 较 好 的 效果 。 然 而 ,这 些 方法 忽 
略 了 SOC 含 量 与 光谱 特征 之 间 的 空间 异 质 性 , 即 假 
定 SOC 含 量 在 不 同 的 空间 位 置 上 对 光谱 特征 的 影 
响 是 不 变 的 。 实 际 上 ,SOC 含量 受 某 一 地 区 土壤 类 
型 .气候 特征 地形 特征 等 自然 因素 和 城市 扩张 . 煤 
VFR .工业 发 展 等 人 为 因素 的 影响 ,使 其 分 布 具 
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有 一 定 的 区 域 性 和 空间 异 质 性 ”。 近 年 来 ,由 于 地 
理 加 权 回 归 (CWR) 能 够 解释 自 变 量 和 因 变 量 之 间 
的 局 部 关系 ,因此 受到 越 来 越 多 的 关注 ”。GCWR 是 
针对 目标 变量 和 解释 变量 之 间 回 归 系 数 的 空间 非 
平稳 性 而 专门 设计 的 , 它 利 用 局 部 数据 局 部 地 测量 
这 些 系数 ”, 对 土壤 属性 空间 分 布 的 预测 已 取得 一 
些 成 效 , 但 将 多 光谱 信息 和 GWR 模型 相 结 合 预测 
SOC 方 面 的 研究 较 少 ,特别 是 关于 西北 旱 区 工业 城 
市 SOC 的 预测 较为 缺乏 。 

石嘴山 市 地 处 西北 干旱 半 干 旱 区 ,生态 环境 极 
其 脆弱 ,而 其 又 是 依托 煤炭 资源 而 发 展 起 来 的 典型 
资源 型 老 工业 城市 。 从 20 世 纪 开 始 , 国 家 就 大 规模 
开发 石嘴山 市 的 煤炭 资源 ,截至 2011 年 石嘴山 市 有 
煤 、 电 钢材、 碳化 硅 、 硅 铁 等 高 耗 能 、 高 污染 企业 
300 多 家 ,三废 排 放量 大 ,因此 影响 了 全 球 气 候 变 
化 和 生态 系统 的 稳定 性 ,使 得 SOC 含量 存在 较 大 的 
空间 变异 性 。 为 此 ,本 文 以 宁夏 回族 自治 区 石嘴山 
市 城市 SOC 为 研究 对 象 , 结 合 GWR 模型 与 遥感 反 
演 方法 ,利用 土壤 样 点 实测 数据 .数字 高 程 模型 
(Digital Elevation Model, DEM) 数 据 和 Landsat 8 多 
光谱 数据 ,通过 提取 地 形 因 子 和 波段 亮度 值 (digital 
number, DN) , 4 Pr SOC 的 空间 变异 规律 ,并 分 别 采 
用 普通 克 里 格 法 (OK) 、 多 元 线性 回归 克 里 格 法 
(RK) 遥感 反 演 方法 (RS) 和 遥感 -地 理 加 权 回 归 克 
里 格 法 (RGCWRK ) 对 研究 区 SOC 含量 进行 预测 , 确 
定 研 究 区 SOC 空 间 分 布 最 优 预 测 方法 ,为 城市 土壤 
属性 预测 提供 数据 支持 。 


1 研究 区 概况 与 数据 来 源 


1.1 研究 区 概况 

石嘴山 市 位 于 宁夏 回族 自治 区 最 北端 ,地 处 宁 
ZR ` 蒙 西 两 个 国家 千 亿 吨 级 煤田 之 间 ,以 黄河 为 界 , 东 
临 鄂尔多斯 台地 , 西 中 银川 平原 北部 '" ,地 理 坐 标 
界 于 106°16' ~ 106?33' E,38?52' ~ 39°05' N。 气 候 
干燥 多 风 , 年 降水 量 180 ~ 190 mm, 2€ f 1.700 ~ 
2 500 mm ,为 典型 的 干旱 半 干 旱 气 候 , 年 平均 气温 
8.4 ~ 9.9 ,年 最 低 平均 气温 -19.4 ~ -23.2 % ,年 最 高 
平均 气温 32.4~ 36.1 CC。 随 着 “一 带 一 路 ”战略 的 提 
出 ,石嘴山 市 成 为 宁夏 内 陆 开放 型 经 济 试验 区 建设 
的 核心 区 之 一 ,但 伴随 着 城市 建设 .工业 .采矿 等 人 为 
因素 的 影响 ,使 石嘴山 市 城市 SOC 表现 出 异 质 性 与 
多 样 性 ,从 而 影响 着 城市 生态 环境 的 稳定 性 。 
1.2 土壤 样本 采集 及 处 理 

2018 年 5 月 在 研究 样 地 内 ,将 研究 区 按 不 同 的 
土地 利用 方式 划分 为 8 大 功能 区 :道路 公园、 工业 
区 .商业 区 .科教 区 .住宅 区 .农田 和 医疗 区 ,通过 手 
持 式 全 球 定位 系统 (GPS) 定 位 ,在 石嘴山 市 范围 内 
按照 不 同 功能 区 随机 选取 151 个 样 点 ,其 中 114 个 
点 用 于 方程 拟 合 ,37 个 点 用 于 方程 验证 。 取 样 深度 
为 0~ 20 cm。 样 本 带 回 实验 室 风 干 .研磨 ,过 100 目 
fifi ETT SOC 含量 测定 ,采用 重 铬 酸 钾 容量 法 一 外 
加 热 法 中 。 人 研究 区 土壤 样 点 分 布 如 图 1 所 示 : 
1.3 影像 数据 及 处 理 

Landsat 8 多 光谱 影像 被 用 作 土 壤 光 谱 信 息 的 


0 6km 


图 1 研究 区 土壤 样 点 分 布 示意 图 


Fig. 1 Distribution of sampling points in study area 
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数据 源 , 按 照 土 壤 裸 露 . 云 量 较 少 、 季 节 相 似 、 影 像 
年 份 邻近 采样 时 间 等 原则 ,选取 了 2018 年 4 月 2 日 
空间 分 辩 率 为 30 m x 30 m 的 多 光谱 影像 。 通 过 
ENVI 5.3 对 影像 进行 辐射 定 标 、 大 气 校正 等 相关 处 
理 , 使 其 最 大 程度 地 与 采样 点 土壤 属性 的 真实 值 接 
近 。 从 遥感 影像 中 所 提取 的 数据 为 DY 值 , 是 传 感 吉 
接收 地 物 反 射 的 能 量 转换 后 的 结果 。 由 于 太阳 高 度 
角 和 大 气 条 件 会 不 可 避免 地 影响 地 面 辐射 ,导致 地 
面 的 电磁 波 波 谱 发 生 畸 变 , 为 了 消除 这 些 因素 的 影 
响 ,需要 对 遥感 图 像 进行 辐射 校正 ,才能 得 到 反射 率 
的 真实 值 2 。 筛 选 出 与 SOC 分 布 相关 的 波段 ,分 析 
表层 SOC 含量 与 1~7 波 段 DV 值 及 DN 值 变换 形式 
的 相关 关系 25 ,建立 SOC 含 量 的 光谱 预测 模型 。 
1.4 地 形 数据 及 处 理 

从 地 理 空间 数据 云 上 下 载 分 辨 率 为 30 m 的 


高 程 /m 
3250 


956 
— d 


DEM 数据 ,利用 ArcGIS 空间 分 析 工 具 , 提取 常 用 的 
地 形 因子 :高 程 坡度 UE In] EHE AES .地 形 起 伏 度 
(图 2) ,其 中 坡度 变 率 是 微分 空间 中 坡度 的 坡度 n] 
面 曲 率 的 信息 在 一 定 程 度 上 可 以 反映 出 来 。 地 形 
起 伏 度 是 指 某 一 特定 区 域 的 最 高 海拔 和 最 低 海拔 
之 间 的 差异 , 它 是 一 个 宏观 数据 ,用 来 描述 区 域 的 
地 形 特 征 。 
1.5 研究 方法 

RK 是 基于 多 元 线性 回归 创建 SOC 与 相关 解释 
变量 之 间 的 关系 ,其 预测 结果 由 预测 模型 计算 得 出 ， 
对 预测 结果 进行 OK 插值 即 为 SOC 的 最 终 预 测 结果 。 

地 理 加 权 回 归 克 里 格 法 (GWRK) 代 替 了 RK 法 
中 的 全 局 拟 合 ,将 其 转换 成 GWR 中 的 局 部 拟 合 , 然 
后 将 插值 点 的 局 部 适用 性 残 差 进行 OK 插值 。 为 了 
更 加 清楚 地 理解 GWRK 法 , 先 对 GWR 模型 做 以 下 


(b) 


起 伏 度 /m 


= 196 
0 


c 湖 


2 研究 区 地 形 特征 图 


Fig.2 Terrain features of the study area 
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解释 : 
Y owr) =B So ) + YB So ) Xi so) 


式 中 : Yu G) 为 因 变 量 了 在 点 s 处 的 模拟 值 ; 
XG O NEk AARETE s 处 的 测量 值 ; so 为 局 
部 估计 的 系数 ,用 加 权 最 小 二 乘法 求 得 ,是 空间 位 
置 的 函数 ,其 权重 矩阵 的 估计 采用 高 斯 函数 ;为 自 
变量 数目 。 

GWRK 法 是 对 CWR 法 的 一 种 深入 推广 , 即 用 
OK 法 对 GCWR 模 型 拟 合 后 的 残 差 进行 插值 ,并 增加 
了 GWR 法 拟 合 的 趋势 ,其 表达 式 为 : 


k 
Y v 59) 7 B4 $9) 十 Yt so) X $5) +e ( So ) 


式 中 : e (30) 为 GWR 模 型 在 点 so 处 拟 合 后 所 得 的 
残 差 ,采用 OK 法 进行 插值 。 

RS 是 基于 实验 室 对 土壤 光谱 反射 特性 的 研究 ， 
在 少量 地 表 数 据 测量 中 ,分 析 遥 感 影像 采样 点 光谱 
特征 的 变化 ,建立 研究 区 土壤 信息 预测 模型 。 大 量 
前 人 研究 表明 SOC 与 光谱 反射 率 存 在 显著 的 负 相 
关 关 系 ' ,利用 这 种 相关 关系 建立 定量 预测 模型 可 
以 较 准 确 模拟 SOC 含 量 。SOC 遥 感 反 演 的 精度 受到 
多 种 因素 影响 ,包括 土壤 其 他 理化 性 质 .地形 因素 
和 大 气 因 素 等 。 首 先 , 对 实测 SOC 与 利用 影像 DN 
值 构建 的 指标 进行 相关 性 分 析 ,得 到 与 SOC 显著 或 
极 显著 相关 的 变量 作为 解释 变量 。 其 次 ,为 避免 数 
据 的 重复 ,假设 SOC 测 量 值 为 因 变量 ,各 波段 DN 值 
为 建 模 样本 的 自 变量 ,利用 SPSS 18.0 软件 ,建立 
SOC 含 量 的 数据 反 演 模型 。 

RGWRK 是 将 多 光谱 遥感 影像 与 GWRK 方 法 相 
结合 ,使 用 遥感 影像 DV 值 与 研究 区 地 形 因 子 来 共 
同 预测 石嘴山 市 城市 SOC 含量 。 

1.6 精度 评价 

随机 生成 的 114 个 用 于 拟 合 的 点 参与 所 有 分 析 
过 程 , 其 余 37 个 测试 点 用 于 验证 空间 搬 值 结果 和 遥 
感 反 演 结 果 。 用 验证 点 估计 值 和 实测 值 的 平均 相 
对 误差 (MRE) 和 均 方 根 误差 (RMSE ) 评 价 不 同方 法 
对 预测 的 准确 性 "。 


SOC 含量 的 实测 值 , n 为 验证 点 个 数 ,在 本 文中 即 
n=37。 当 MRE 和 RMSE 接 近 于 0 时 ,插值 精度 较 高 。 
以 RMSE 为 指标 ,比较 不 同 插值 方法 与 参考 方 
法 的 相对 精度 提高 值 (RD)" ,公式 如 下 : 
pr- RUSEn RMSE; 
RMSE, 
式 中 : RMSE, 为 参照 方法 的 RMSE , RMSE, 为 其 余 
方法 的 RMSE 。 
1.7 数据 处 理 
利用 SPSS 18.0 软 件 对 样本 数据 进行 描述 性 统 
计 和 多 元 线性 回归 ,并 且 进 行 Kolmogorov-Smirnov 
(K-S) 检 验 ,对 数据 进行 OK 插值 的 前 提 条 件 是 其 要 
符合 正 态 分 布 ""。 在 GS+7.0 软件 中 完成 数据 的 半 
方差 分 析 , 在 GWR 4.0 中 完成 GWR 分 析 , 用 ArcGIS 
10.3 软件 分 别 对 研究 区 的 SOC 进行 OK 插值 .RK 插 
值 和 RGWRK 插 值 。 


2 结果 与 分 析 


2.1 土壤 有 机 碳 描述 性 统计 

由 表 1 可 知 , 不 同 功能 区 SOC 含量 大 小 表现 为 : 
工业 区 (23.40 g* kg )> 医 疗 区 (21.18 g* kg)» T MV DX. 
(20.02 g- kg )> 道 路 (19.29 g:kg')> 住 宅 区 (16.70 g- 
kg )»73 pä (13.87 g-kg 和)> 农 田 (12.16 gkg”) >F} 
[X (11.14 g* kg) ,其 中 工业 区 是 农田 的 1.92 倍 。 不 
同 功 能 区 SOC 含 量 与 乌鲁木齐 市 “研究 结果 接近 ， 
但 与 开封 市 ”、 南 京 市 有 一 定 的 差异 ;由 于 各 功 
能 区 的 划分 以 及 土壤 的 采样 密度 、 计 算 单位 有 差 
异 ,使 得 不 同城 市 不 同 功 能 区 SOC 空间 分 布 特征 不 
确定 。 变 异 系 数 的 大 小 表示 土壤 属性 空间 变异 性 
的 大 小 ,根据 CV 评价 准则 ,CV 小 于 10% 时 变异 性 
较 弱 ,CV 大 于 100% 时 变异 性 较 强 , 位 于 二 者 之 间 的 
属于 中 等 变异 ,本 研究 区 各 功能 区 SOC 的 CV 在 
33.90% ~76.75% 之 间 变 化 , 均 属 于 中 等 程度 变异 ， 
说 明 在 不 同 功能 区 SOC 分 布 不 均匀 。 各 功能 区 SOC 
含量 均 通 过 K-S 检 验 呈 (p > 0.05) ,符合 正 态 分 布 。 

对 8 个 功能 区 SOC 含 量 进行 方差 分 析 ,结果 表 
明 不 同 功能 区 SOC 存在 显著 差异 (p < 0.05)。 工 业 
区 与 公园 .科教 区 和 农田 存在 显著 差异 (p < 0.05 ) 。 
科教 区 与 道路 存在 显著 差异 (p < 0.05) 。 医 疗 区 与 
科教 区 和 农田 显著 差异 (p < 0.05)。 说 明 土 地 利用 
格局 是 影响 石嘴山 市 城市 SOC 空间 分 布 的 重要 因 
素 ,也 说 明 研 究 区 SOC 易 受 交 通 ,城市 建设 和 商业 


x 100% 
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表 1 SOC 描述 统计 量 
Tab.1 Descriptive statistics of SOC 


功能 区 样 点 数 最 小 值 最 大 值 土壤 有 机 碳 偏 度 峰 度 变异 系数 /% K-S 
Min. Max. SOC Skewness Kurtosis CV p (& 
道路 16 6.17 26.00 19.29+6.54abc -0.55 -0.92 33.90 0.453 
工业 区 25 5.18 66.92 23.40+17.96a 1.32 0.73 76.75 0.065 
商业 区 6 4.00 26.29 20.02+8.47abcd -1.78 3.12 42.31 0.789 
医疗 区 12 3.65 33.60 21.18+9.78ab -0.35 -0.78 46.18 0.994 
住宅 区 15 3.73 34.24 16.70x9.43abcd 0.29 -0.72 56.47 0.996 
公园 10 5.60 25.87 13.87+6.11bcd 0.67 0.56 44.05 0.764 
科教 区 14 1.31 26.78 11.14+8.01d 0.922 -0.05 71.90 0.856 
农田 16 3.09 37.31 12.16+7.98cd 2.23 6.49 65.63 0.394 
合计 114 1.31 66.92 17.61x11.67 1.67 4.28 66.27 0.068 


活动 等 人 类 活动 影响 , 且 在 城市 空间 分 布 上 存在 较 
大 差异 。 
22 土壤 有 机 碳 空 间 变 异 特征 分 析 

利用 地 统计 学 半 方 差 理论 分 析 SOC 含量 的 空 
间 变 异 结构 。 表 2 显示 了 4 种 模型 对 SOC 半 变异 函 
数 的 模拟 结果 及 交叉 检验 结果 ,根据 各 模型 交 又 验 
证 结果 ,高 斯 模型 的 相关 系数 平方 和 (及 = 0.694) 最 
大 , 且 回 归 系 数 标准 误差 (Se = 0.687) 和 标准 误差 预 
ill (& (Sep = 0.097) 在 四 个 模型 中 最 小 。 因 此 ,该 
SOC 含 量 的 最 佳 拟 合 模型 为 高 斯 模型 。 高 斯 模型 中 
基 台 值 (Co+C) 为 116.10, 表 示 变 量 的 最 大 变异 程 
度 。GoV(CotC) 反 映 了 SOC 的 空间 依赖 性 3。 一 般 
来 说 ,比值 <0.25 时 空间 变异 性 表现 为 强 空 间 相 关 
性 , 介 于 0.25 ~ 0.75 时 表现 为 中 等 强度 空间 相关 性 ， 
比值 >0.75 时 表现 为 弱 空 间 相 关 性 。 变 程 表明 了 人 研 
究 样 本 区 域 的 相似 程度 ,可 作为 确定 环境 因素 之 间 
距离 的 一 种 度量 六 , 变 程 越 小 ,表明 SOC 在 空间 上 
的 相互 作用 距离 越 小 ;反之 , 变 程 越 大 ,变量 相互 依 
存 的 空间 距离 就 越 大 , 即 变异 性 越 小 ,本 研究 SOC 
的 变 程 为 2 424.87 mo 


表 2 SOCZE 


注 :土壤 有 机 碳 数 据 为 平均 值 + 平 均 误 差 , 不 同 功能 区 不 同 小 写字 母 表 示 差 异 显著 (p < 0.05) 


2.3 影像 DN 值 \ 地 形 因 子 与 土壤 有 机 碳 含量 相关 性 

分 析 

对 研究 区 SOC 分别 与 遥感 图 像 各 波段 DN 值 .地 
形 指标 进行 Pearson 相关 分 析 ( 表 3)。 由 表 3 可 以 看 
出 ,B1-B7 与 B2-B7 之 间 DN 值 的 相关 系数 最 高 (r = 
0.993,p < 0.01); 其 次 为 B2-B6 与 B3-B6(r = 0.979, 
p< 0.01) ,表明 这 些 波 段 DN 值 之 间 的 重 羡 较 多 。 
SOC 与 各 波段 DN 值 之 差 存 在 显著 相关 关系 (p < 
0.01) ,其 中 1、2、3、4、5 波段 与 6.7 波 段 之 差 与 SOC 
表现 出 极 显 著 负 相关 (p < 0.01) ,相关 系数 达到 0.260 
及 以 上 ,其 中 B4-B7 与 SOC 相关 系数 最 高 。 同 时 ， 
SOC 与 高 程 呈 极 显著 正 相 关 关 系 (r = 0.321,p < 
0.01) ,表明 海拔 越 高 ,SOC 含量 越 高 。SOC 与 坡度 也 
旦 极 显 著 正 相关 关系 (r = 0.346,p < 0.01) ,说 明 在 一 
定 范 围 内 ,坡度 越 大 处 ,有 机 碳 越 容易 富 集 。 此 外 ， 
SOC 与 地 表 起 伏 度 的 相关 性 也 达到 了 显著 正 相关 水 
平 r= 0.312,p < 0.01) ,反映 了 地 形 也 会 对 SOC 含量 
产生 影响 。 
2.4 土壤 有 机 碳 多 元 线性 回归 预测 模型 

逐步 回归 不 仅 可 以 保证 与 SOC 显著 相关 的 自 


异 模型 参数 


Tab. 2 Semivariance parameters of SOC 


标准 误差 预测 值 


模型 块 金 值 C[。 基 台 值 CotC 变 程 4 Co / (CC) DOES — ”标准 误差 Se 
Ps ep 
线性 模型 95.08 123.82 11217.45 0.77 0.247 0.723 0.124 
sot 球状 模型 43.00 116.00 2850.00 0.37 0.690 0.705 0.125 
` 指数 模型 54.70 116.00 2220.00 0.47 0.658 0.709 0.129 
高 斯 模型 2.10 116.10 2424.87 0.02 0.694 0.687 0.097 
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表 3 SOC 与 DN 值 差 值 \ 地 形 因子 的 相关 性 分 析 结 果 


Tab.3 Correlation between SOC and the difference of DN and terrain factor 


SOC B1-B7 B2-B6 B2-B7 B3-B6 


B3-B7 B4-B6 B4-B7  B5-B6 高 程 坡度 


SOC 1 

B1-B7 -0.251™ 1 

B2-B6 -0.279" 0.967” 1 

B2-B1 -0.281™ 0.993™ 0.973 1 

B3-B6 -0.266* 0.957” 0.9797 0.957” 1 
B3-B7 -0.2607 0.966” 0.9247 0.966" 0.960” 
B4-B6 -0.343" 0.890” 0.936” 0.912” 0.964” 
B4-B7 -0.349" 0.865" 0.840” 0.890” 0.887” 
B5-B6 -0.312^ 0.781” 0.727” 0.794” 0.763” 


坡度 0.346” -0.145 -0.170  -0.154 -0.163  -0.139  -0.182 -0.154 
起 伏 度 0.312” -0.139  -0.171  -0.147  -0.165  -0.132  -0.184 -0.147 


变量 进入 回归 模型 ,而 且 可 以 去 除 自 变量 之 间 的 共 
线性 下 。 本 人 研究 以 还 感 影像 各 波段 DN 值 以 及 地 形 
因子 为 自 变量 ,对 研究 区 SOC 进行 多 元 线性 逐步 回 
归 分 析 , 模 型 1 如 下 : 
SOC=-0.018x(B1-B7)+0.047x(B3-B7)+ 
1.081 x Slope — 0.040 x (B4 — B7) + 63.946 
式 中 ;SOC 为 土壤 有 机 碳 含 量 ;B1-B7 为 波段 1 与 波 
Et7 DN (E25 (8 ; B3-B7 为 波段 3 与 波段 7 DN 值 差 
值 ;B4-B7 为 波段 4 与 波段 7 DN (8 25 TH ; Slope 为 研 
由 模型 可 知 ,波段 1、3、4 和 波段 7 的 差 值 与 地 
乡 因子 坡度 是 研究 区 SOC 多 元 线性 回归 分 析 的 最 
佳 辅助 变量 。 回 归 模 型 的 决定 系数 尼 相 对 较 低 , 仅 
为 0.275。 之 前 的 研究 中 也 多 次 出 现 过 类 似 的 结 
如 张 素 梅 等 “在 预测 农安 县 SOM 分布 时 ,回归 方程 
的 决定 系数 为 0.220, 刘 焕 军 等 在 预测 东北 典型 黑 
土 区 田 块 土壤 有 机 质 进 行 逐 步 回 归 拟 合 时 ,回归 方 
程 的 决定 系数 为 0.341。 在 ArcGIS 10.3 中 对 通过 多 
元 线性 回归 模型 得 到 的 预测 值 进行 OK 插值 ,得 到 
最 终 预测 结果 ,并 进行 精度 验证 。 
2.5 土壤 有 机 碳 遥 感 反 演 结果 分 析 
利用 SPSS 软件 对 SOC 含 量 与 各 波段 DN 值 进行 
多 元 线性 逐步 回归 分 析 ,构建 得 到 模型 2 如 下 : 
SOC = -0.021 x (B1 - B7) + 0.052 x (B3 - B7) - 
0.043 x (B4 - B7) + 73.836 
式 中 ;SOC 为 土壤 有 机 碳 含 量 ;B1-B7 为 波段 1 与 波 


(1) 


(2) 


0.925™ 1 

0.941” 0.9942" 1 
0.835" 0.760” 0.827” 1 
高 程 0.3211 -0.114 -0.136 -0.125 -0.138 20.123 -0.184 -0.175 


-0295" 1 
-0.257' 0.736" 1 
-0.241” 0.747" 0.985" 1 


注 : B; 表 示 Landsat 8 影像 第 i 波段 DN 值 ,BB 表示 第 i 与 第 j 波 段 DN 值 之 差 ;“ 表 示 在 0.01 水 平 上 显著 相关 ,表示 在 0.05 水 平 上 显著 相关 


Et7 DN (82518 ; B3- B7 为 波段 3 与 波段 7 DN 值 差 
值 ;B4-B7 为 波段 4 与 波段 7DN 值 差 值 。 

根据 模型 可 知 ,模型 以 (B1-B7)、(B3-B7)、(B4- 
B7) 作 为 自 变 量 , 反 演 模 型 拟 合 效果 较 好 。 逐 步 回归 
方法 可 以 消除 数据 之 间 的 宛 余 , 对 于 SO0C 含量 的 反 
演 而 言 ,用 于 建 模 各 DN 值 差 值 之 间 存 在 着 重复 信 
息 。 回 归 方 程 的 系数 检验 结果 见 表 4, 波 段 4 与 波段 
7 的 DN 值 差 值 系数 通过 了 统计 学 显著 性 检验 (p < 
0.001) ,但 波段 1,3 与 波段 7DN 值 差 值 的 系数 未 通 
过 检验 ,p 值 (0.003 和 0.029) 较 小 ,因此 ,在 方程 中 仍 
保留 了 (B83-B7) 和 (B1-B7)。 


表 4 回归 方程 系数 的 检验 


Tab.4 Test for coefficients of the regression equation 


系数 系数 值 ”标准 误差 1 检验 概率 
常数 73.836 12.917 5.716 P<0.001 
B4-B7 -0.043 0.010 -4.261 P«0.001 
B3-B7 0.052 0.017 3.042 P-0.003 
B1-B7 -0.021 0.009 -2.216 P=0.029 


2.6 土壤 有 机 碳 遥 感 -地 理 加 权 回 归 克 里 格 插值 
GWR 是 针对 目标 变量 和 解释 变量 之 间 回 归 系 
数 的 空间 非 平 稳 性 而 专门 设计 的 , 它 利 用 局 部 数据 
局 部 测量 这 些 系 数 。 与 传统 回归 方法 不 同 ,GWR 法 
在 空间 上 不 是 全 局 的 (无 变化 ) ,而 是 依赖 于 不 同 空 
间 位 置 的 。GWRK 是 一 个 混合 技术 ,并且 是 CWR 
方法 的 延续 。GWRK 是 由 确定 和 随机 两 个 部 分 组 
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成 ,这 两 个 部 分 分 别 建 模 。 确 定性 的 组 成 部 分 是 由 
SOC 来 建 模 并 且 使 用 可 用 的 相关 变量 信息 来 预测 目 
标 变量 的 趋势 。 随 机 的 组 成 部 分 ( 残 差 ) 可 以 用 OK 
插值 ,并 添加 到 估计 的 趋势 。 残 差 被 认为 可 能 具有 
一 些 可 以 建 模 的 空间 相关 结构 , 残 差 是 模型 的 组 成 
部 分 ,不 能 够 被 确定 性 成 分 解释 ,但 重要 的 是 要 补 
充 , 因 为 它 能 够 帮助 解释 目标 变量 在 空间 上 的 变 
化 。GWRK 方 法 能 够 用 于 估计 残 差 ,从 而 更 有 效 地 
预测 趋势 。 为 了 方便 比较 ,应 用 GWRK 模 型 时 , 选 
取 相 同 的 波段 信息 (B1-B87) ,(B3- B7) 、(B4-B7) 和 
地 形 因 子 高 程 .坡度 .起 伏 度 进行 建 模 , 通 过 使 用 
GWR 4.0 软 件 来 构建 GWR 模型 ,并 通过 MRE、RMSE 
和 RI 对 其 进行 精度 验证 。 
2.7 4 种 不 同 预测 方法 下 土壤 有 机 碳 空 间 分 布 特征 
采用 四 种 方法 分 别 对 SOC 含量 的 空间 分 布 进 
行 OK 插 值 ,得 出 SOC 含 量 的 空间 分 布 图 。 从 分 布 
特征 看 ,4 种 方法 预测 的 SOC 含量 空间 分 布 具有 相 
似 的 格局 , 即 研究 区 北部 SOC 含量 要 大 于 南部 ,高 


SOC / g'kg! 
m 7.44~10.03 


(b) RK aa! 


0 mm 10.04-12.62 
m 12.63-15.21 


m 25.58-28.16 
m 28.17-30.75 


SOC / gkg' 

m 8.46~10.22 
æ 10.23-11.98 
zu 11.99-13.74 
= 13.75-15.49 
= 15.5-17.25 
= 17.26-19.01 
æ 19.02-20.76 
um 20.77-22.52 
m 22.53-24.28 


mE 


值 区 分 布 在 正 北 和 东南 部 。 从 制图 效果 看 ,OK 法 
得 到 的 图 较为 粗略 , 斑 块 非常 零碎 ,不 符合 现实 中 
SOC 含量 空间 渐变 的 特征 ;RK 法 和 RS 法 的 插值 结 
果 体现 了 SOC 含 量 随 着 各 变量 的 变化 情况 ,但 是 也 
有 人 少 部 分 斑 块 的 边界 较为 零碎 , 且 RS 法 预测 出 的 
SOC 含量 比 起 其 他 方法 偏 低 ;RGCWRK 法 的 SOC 含 
量 预测 图 在 空间 变化 上 表现 为 自然 .平缓 的 过 渡 ， 
有 效 避 免 了 OK 法 和 RK 法 图 斑 边界 两 侧 的 破碎 程 
ER 突变 明显 的 缺陷 ,还 反映 出 SOC 含 量 的 空间 
渐变 特征 , 斑 块 边界 更 为 光滑 ,感官 上 更 符合 SOC 
含量 的 实际 分 布 状 况 。 
2.8 土壤 有 机 碳 预测 方法 精度 比较 

通过 计算 37 个 验证 点 的 MRE、RMSE 和 RI, 检 
验 不 同方 法 的 预测 精度 。 从 表 5 中 可 知 , 引 入 地 形 
子 和 波段 DN 值 差 值 的 RGWRK 法 各 项 精度 指标 
均 优 于 OK 法 、RK 法 和 RS 法 。 分 别 以 OK 法 、RK 法 
和 RS 法 作为 参照 ,在 预测 SOC 的 精度 上 分 别 相对 提 
高 了 10.05% ,8.7996 48.9296 , t RGWRK 法 是 最 适 


N 
人 SOC/gkg 


9 


m 7.01~9.87 
um 9.88~12.73 
zu 12.74~15.59 
r3 15.6~18.44 
c3 18.45~21.3 
r3 21.31-24.16 
m 24.17-27.02 
um 27.03-29.88 
m 29.89-32.74 
= 


s 3b 


SOC / g kg! 
— . m 16.19-20.14 
à m 20.15-24.09 
pe = 24.1~28.04 
” r328.05-31.99 
r332-35.94 
= 35.95-39.89 
m 39.9—43.84 
m 43.85-47.79 
m 47.8—51.74 


= 


图 3 4 种 方法 空间 插值 结果 图 


Fig. 3 Spatial distribution of SOC by four methods 


表 S$ SOC 的 不 同 预测 方法 精度 比较 


Tab. 5 Accuracy assessment of predicted SOC contents by using different methods 


预测 方法 平均 相对 误差 均 方 根 误差 相对 提高 度 /% 相对 提高 度 /% 相对 提高 度 /% 
MRE RMSE RI RI RI 

OK 0.517 19.694 - - - 

RK 0.505 19.423 1.38 - 0.14 

RS 0.507 19.451 1.23 - - 

RGWRK 0.412 17.715 10.05 8.79 8.92 
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合 本 研究 区 的 最 优 插值 方法 。 
3 讨论 


通过 对 石嘴山 市 城市 不 同 功 能 区 SOC 含量 特 
征 分 析 显 示 ,城市 SOC 含量 与 人 为 活动 密切 相关 ， 
石嘴山 市 表层 SOC 平 均 含 量 为 17.61 g'kg ,属于 中 
等 水 平 ,与 开封 市 "(19.10 g* kg) Bp atri (16.29 
gkg ) 等 城市 含量 相近 , 略 高 于 上 海 市 *(11.91 ge 
kg) .杭州 市 (13.94 g*kg') 等 东部 城市 ,明显 低 于 
乌鲁木齐 市 "(33.05 gkg”)。 罗 上 华 33 对 北京 市 绿 
地 土壤 碳 氮 的 研究 表明 ,城市 化 与 工业 化 加 快 了 城 
市 土壤 碳 库 的 累积 。 石 嘴 山 市 现 处 于 城市 化 发 展 
时 期 ,人口 数 量 不 断 增 长 ,道路 .城市 建筑 不 断 施 
工 ,使 得 石嘴山 市 城市 SOC 含量 有 不 同 程度 的 富 
集 。 从 分 布 特征 看 ,研究 区 北部 的 SOC 含量 要 大 于 
南部 ,结合 采样 点 的 功能 区 分 区 可 知 北部 和 东南 部 
是 石嘴山 市 工业 园区 ,工业 区 的 SOC 含量 较 高 ,而 
南部 是 农田 ,其 SOC 含 量 较 低 。 研 究 区 SOC 呈现 不 
规则 分 布 与 其 土地 利用 方式 和 人 类 活动 有 着 密切 
关系 ,其 中 工厂 排放 的 “三 废物 质 成 为 影响 石嘴山 
Ti SOC 的 主要 原因 之 一 ,同时 影响 着 城市 碳 循环 与 
生态 系统 的 稳定 。 与 乌鲁木齐 市 北京 市 ”、 南 京 
市 中 开封 市 和 小、 上 海 市 中、 杭州 市 池 研 究 结 果 相 
一 致 。 

OK 法 是 目前 土壤 属性 空间 插值 最 具 代 表 性 的 
方法 。 该 方法 利用 采样 点 的 结构 特点 ,给 出 了 插值 
结果 的 佑 计 误差 ,但 没有 考虑 土壤 属性 与 环境 因素 
的 密切 联系 , 且 插 值 精度 取决 于 采样 密度 等 因素 ， 
因此 在 环境 复杂 的 地 区 适应 性 较 差 。 基 于 权 函 数 
灵活 变化 的 GWR 模型 ,克服 了 回归 模型 中 不 同 空 
间 位 置 的 环境 因子 权重 一 致 的 不 足 ,实现 了 局 部 空 
间 信 息 的 捕获 ,并 将 OK 法 对 残 差 信息 的 充分 利用 
相 结 合 , 提 高 了 插值 精度 。MISHRA 45 AP 34 GWR 
与 RK 法 进行 比较 ,利用 一 些 地 形 属性 、 气 候 数 据 、 
土地 利用 数据 和 NDVI 数 据 作为 解释 变量 ,预测 美 
国 中 西部 7 个 州 的 SOC 空间 分 布 ,他 们 发 现 ,GWR 
在 捕捉 变异 性 和 提高 准确 性 方面 略 优 于 RK 法 。 而 
众多 学 者 研究 发 现 SOC 含量 和 土壤 反射 率 高 低 有 
明显 关系 3”, 本 文 利用 RGWRK 法 预测 的 石嘴山 城 
市 SOC 含 量 的 精度 最 高 ,最 符合 实际 情况 。 因 此 , 利 
用 多 光谱 遥感 影像 与 GWR 模型 相 结合 的 方法 预测 


城市 SOC 是 一 次 有 益 尝 试 。 同 时 影响 城市 SOC 含量 
预测 结果 准确 性 的 因素 很 多 ,如 区 域 差异 .土壤 水 
分 .气候 因素 以 及 遥感 影像 的 空间 分 辩 率 等 。 后 续 
人 研究 应 当 从 多 方面 考虑 影响 SOC 变化 的 因素 ,寻找 
相关 的 土壤 光谱 特征 参数 ,如 反射 率 的 数学 变换 形 
式 .波段 组 合 等 ,可 以 提高 遥感 监测 SOC 变化 的 
精度 。 


4 结论 


vd 


以 宁夏 石嘴山 市 各 功能 区 为 研究 区 ,以 裸 土 期 
Landsat 8 影像 和 30 m YEK DEM 影像 为 基础 数 
据 ,利用 地 统计 方法 和 RS 法 对 研究 区 SOC 进行 预 
测 ,得 出 以 下 结论 : 

(1) 石 嘴 山 市 SOC 空间 变异 系数 为 66.27% , Œ 
中 等 变异 水 平 ,SOC 含 量 范围 为 1.31 ~ 66.92 grkg”, 
变化 较 广 ,不 同 功能 区 SOC 含量 差异 显著 ,表现 为 ; 
工业 区 > 医疗 区 > 商业 区 > 道路 > 住宅 区 > 公园 > 农田 > 
科教 区 。 

(2) 根 据 地 形 因子 与 SOC 相关 性 分 析 可 知 ,地 
形变 化 和 遥感 影像 波段 信息 同时 影响 着 研究 区 SOC 
含量 的 空间 分 布 。 相 较 于 纯粹 的 遥感 波段 信息 ,对 
波段 DN 值 进行 数学 变换 得 到 的 数据 与 SOC 相关 性 
更 强 。 

(3) 通 过 对 4 种 方法 的 精度 进行 比较 ,引入 地 形 
因子 和 波段 DN 值 差 值 的 RGWRK 法 MRE 和 RMSE 
均 为 最 低 。 分 别 以 OK 法 .RK 法 和 RS 法 作为 参照 ， 
在 预测 SOC 的 精度 上 分 别 相 对 提高 了 10.05% 、 
8.79% 和 8.92%, 故 RGWRK 法 是 最 适合 本 研究 区 的 
最 优 插值 方法 。 
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Spatial distribution of soil organic carbon in Shizuishan based on 
multispectral and geographically weighted regression model 


XIAZi-shu", BAIYi-cu", BAO Wei-bin", ZHONG Yan-xia, WANG You-qi" 
(1 College of Resources and Environment, Ningxia University, Yinchuan 750021, Ningxia, China; 2 Arid Area Characteristic 


Resources and Environmental Governance Department of Education International 


Cooperation Joint Laboratory, Yinchuan 150021, Ningxia, China) 


Abstract: The distribution of soil organic carbon (SOC) in urban areas is influenced by human factors such as ur- 
ban construction and industrialization, and thus shows obvious spatial differences. To reveal the impact of human ac- 
tivities such as urbanization and industrialization on SOC in Shizuishan City, the spatial distribution of SOC in Shi- 
zuishan was predicted using ordinary Kriging (OK), multiple linear regression Kriging (RK), remote sensing inver- 
sion (RS), and remote sensing- geographically weighted regression Kriging (RGWRK). The SOC content changed 
from 1.31 g/kg to 66.92 g* kg ', with an average of 17.61 g* kg '. There were significant differences in SOC content 
among different functional areas in Shizuishan (p « 0.05). Specifically, the performance in decreasing order was, in- 
dustrial areas, medical areas, commercial areas, roads, residential areas, parks, farmlands, scientific and education- 
al areas. The coefficient of variation of SOC content was 66.2796 and there was moderate variation. The best fitting 
model was the Gauss model, and Cy/(C,-* C) was 0.02, indicatingstrong spatial autocorrelation. There was a signifi- 
cant correlation (p < 0.01) between SOC and differences of DN values in different bands of RS images (B1-B7, B3- 
B7, and B4— B7) and topographic factors (elevation, slope, and relief amplitude). Comparing the results of the four 
methods demonstrated that SOC predicted by RGWRK was the most accurate, being 10.0596 higher than that of 
OK, 8.7996 higher than that of RK, and 8.92% higher than that of RS. SOC content in the northern part of Shizuis- 
han was higher than the southern part. The SOC content in industrial areas was 1.92 times higher than farmlands, in- 
dicating that the SOC content in the urban area was enriched. 

Key words: soil organic carbon; spatial prediction; remote sensing inversion; geographically weighted regression 
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